3.1.97 \(\int x (A+B x^2) (a+b x^2+c x^4)^3 \, dx\) [97]

3.1.97.1 Optimal result
3.1.97.2 Mathematica [A] (verified)
3.1.97.3 Rubi [A] (verified)
3.1.97.4 Maple [A] (verified)
3.1.97.5 Fricas [A] (verification not implemented)
3.1.97.6 Sympy [A] (verification not implemented)
3.1.97.7 Maxima [A] (verification not implemented)
3.1.97.8 Giac [A] (verification not implemented)
3.1.97.9 Mupad [B] (verification not implemented)

3.1.97.1 Optimal result

Integrand size = 23, antiderivative size = 166 \[ \int x \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{2} a^3 A x^2+\frac {1}{4} a^2 (3 A b+a B) x^4+\frac {1}{2} a \left (a b B+A \left (b^2+a c\right )\right ) x^6+\frac {1}{8} \left (3 a B \left (b^2+a c\right )+A \left (b^3+6 a b c\right )\right ) x^8+\frac {1}{10} \left (b^3 B+3 A b^2 c+6 a b B c+3 a A c^2\right ) x^{10}+\frac {1}{4} c \left (b^2 B+A b c+a B c\right ) x^{12}+\frac {1}{14} c^2 (3 b B+A c) x^{14}+\frac {1}{16} B c^3 x^{16} \]

output
1/2*a^3*A*x^2+1/4*a^2*(3*A*b+B*a)*x^4+1/2*a*(a*b*B+A*(a*c+b^2))*x^6+1/8*(3 
*a*B*(a*c+b^2)+A*(6*a*b*c+b^3))*x^8+1/10*(3*A*a*c^2+3*A*b^2*c+6*B*a*b*c+B* 
b^3)*x^10+1/4*c*(A*b*c+B*a*c+B*b^2)*x^12+1/14*c^2*(A*c+3*B*b)*x^14+1/16*B* 
c^3*x^16
 
3.1.97.2 Mathematica [A] (verified)

Time = 0.05 (sec) , antiderivative size = 154, normalized size of antiderivative = 0.93 \[ \int x \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{560} x^2 \left (280 a^3 A+140 a^2 (3 A b+a B) x^2+280 a \left (a b B+A \left (b^2+a c\right )\right ) x^4+70 \left (3 a B \left (b^2+a c\right )+A \left (b^3+6 a b c\right )\right ) x^6+56 \left (b^3 B+3 A b^2 c+6 a b B c+3 a A c^2\right ) x^8+140 c \left (b^2 B+A b c+a B c\right ) x^{10}+40 c^2 (3 b B+A c) x^{12}+35 B c^3 x^{14}\right ) \]

input
Integrate[x*(A + B*x^2)*(a + b*x^2 + c*x^4)^3,x]
 
output
(x^2*(280*a^3*A + 140*a^2*(3*A*b + a*B)*x^2 + 280*a*(a*b*B + A*(b^2 + a*c) 
)*x^4 + 70*(3*a*B*(b^2 + a*c) + A*(b^3 + 6*a*b*c))*x^6 + 56*(b^3*B + 3*A*b 
^2*c + 6*a*b*B*c + 3*a*A*c^2)*x^8 + 140*c*(b^2*B + A*b*c + a*B*c)*x^10 + 4 
0*c^2*(3*b*B + A*c)*x^12 + 35*B*c^3*x^14))/560
 
3.1.97.3 Rubi [A] (verified)

Time = 0.41 (sec) , antiderivative size = 164, normalized size of antiderivative = 0.99, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.130, Rules used = {1576, 1140, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int x \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx\)

\(\Big \downarrow \) 1576

\(\displaystyle \frac {1}{2} \int \left (B x^2+A\right ) \left (c x^4+b x^2+a\right )^3dx^2\)

\(\Big \downarrow \) 1140

\(\displaystyle \frac {1}{2} \int \left (B c^3 x^{14}+c^2 (3 b B+A c) x^{12}+3 c \left (B b^2+A c b+a B c\right ) x^{10}+\left (B b^3+3 A c b^2+6 a B c b+3 a A c^2\right ) x^8+\left (3 a B \left (b^2+a c\right )+A \left (b^3+6 a c b\right )\right ) x^6+3 a \left (a b B+A \left (b^2+a c\right )\right ) x^4+a^2 (3 A b+a B) x^2+a^3 A\right )dx^2\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {1}{2} \left (a^3 A x^2+\frac {1}{2} a^2 x^4 (a B+3 A b)+\frac {1}{2} c x^{12} \left (a B c+A b c+b^2 B\right )+a x^6 \left (A \left (a c+b^2\right )+a b B\right )+\frac {1}{5} x^{10} \left (3 a A c^2+6 a b B c+3 A b^2 c+b^3 B\right )+\frac {1}{4} x^8 \left (A \left (6 a b c+b^3\right )+3 a B \left (a c+b^2\right )\right )+\frac {1}{7} c^2 x^{14} (A c+3 b B)+\frac {1}{8} B c^3 x^{16}\right )\)

input
Int[x*(A + B*x^2)*(a + b*x^2 + c*x^4)^3,x]
 
output
(a^3*A*x^2 + (a^2*(3*A*b + a*B)*x^4)/2 + a*(a*b*B + A*(b^2 + a*c))*x^6 + ( 
(3*a*B*(b^2 + a*c) + A*(b^3 + 6*a*b*c))*x^8)/4 + ((b^3*B + 3*A*b^2*c + 6*a 
*b*B*c + 3*a*A*c^2)*x^10)/5 + (c*(b^2*B + A*b*c + a*B*c)*x^12)/2 + (c^2*(3 
*b*B + A*c)*x^14)/7 + (B*c^3*x^16)/8)/2
 

3.1.97.3.1 Defintions of rubi rules used

rule 1140
Int[((d_.) + (e_.)*(x_))^(m_.)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x 
_Symbol] :> Int[ExpandIntegrand[(d + e*x)^m*(a + b*x + c*x^2)^p, x], x] /; 
FreeQ[{a, b, c, d, e, m}, x] && IGtQ[p, 0]
 

rule 1576
Int[(x_)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^( 
p_.), x_Symbol] :> Simp[1/2   Subst[Int[(d + e*x)^q*(a + b*x + c*x^2)^p, x] 
, x, x^2], x] /; FreeQ[{a, b, c, d, e, p, q}, x]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.1.97.4 Maple [A] (verified)

Time = 0.13 (sec) , antiderivative size = 170, normalized size of antiderivative = 1.02

method result size
norman \(\frac {a^{3} A \,x^{2}}{2}+\left (\frac {3}{4} A \,a^{2} b +\frac {1}{4} B \,a^{3}\right ) x^{4}+\left (\frac {1}{2} A c \,a^{2}+\frac {1}{2} A a \,b^{2}+\frac {1}{2} B \,a^{2} b \right ) x^{6}+\left (\frac {3}{4} A a b c +\frac {1}{8} A \,b^{3}+\frac {3}{8} a^{2} B c +\frac {3}{8} B a \,b^{2}\right ) x^{8}+\left (\frac {3}{10} A a \,c^{2}+\frac {3}{10} A \,b^{2} c +\frac {3}{5} B a b c +\frac {1}{10} B \,b^{3}\right ) x^{10}+\left (\frac {1}{4} A b \,c^{2}+\frac {1}{4} B a \,c^{2}+\frac {1}{4} B \,b^{2} c \right ) x^{12}+\left (\frac {1}{14} A \,c^{3}+\frac {3}{14} B b \,c^{2}\right ) x^{14}+\frac {B \,c^{3} x^{16}}{16}\) \(170\)
gosper \(\frac {1}{2} a^{3} A \,x^{2}+\frac {3}{4} x^{4} A \,a^{2} b +\frac {1}{4} x^{4} B \,a^{3}+\frac {1}{2} x^{6} A c \,a^{2}+\frac {1}{2} x^{6} A a \,b^{2}+\frac {1}{2} x^{6} B \,a^{2} b +\frac {3}{4} x^{8} A a b c +\frac {1}{8} x^{8} A \,b^{3}+\frac {3}{8} x^{8} a^{2} B c +\frac {3}{8} x^{8} B a \,b^{2}+\frac {3}{10} x^{10} A a \,c^{2}+\frac {3}{10} x^{10} A \,b^{2} c +\frac {3}{5} x^{10} B a b c +\frac {1}{10} x^{10} B \,b^{3}+\frac {1}{4} x^{12} A b \,c^{2}+\frac {1}{4} x^{12} B a \,c^{2}+\frac {1}{4} x^{12} B \,b^{2} c +\frac {1}{14} x^{14} A \,c^{3}+\frac {3}{14} x^{14} B b \,c^{2}+\frac {1}{16} B \,c^{3} x^{16}\) \(194\)
risch \(\frac {1}{2} a^{3} A \,x^{2}+\frac {3}{4} x^{4} A \,a^{2} b +\frac {1}{4} x^{4} B \,a^{3}+\frac {1}{2} x^{6} A c \,a^{2}+\frac {1}{2} x^{6} A a \,b^{2}+\frac {1}{2} x^{6} B \,a^{2} b +\frac {3}{4} x^{8} A a b c +\frac {1}{8} x^{8} A \,b^{3}+\frac {3}{8} x^{8} a^{2} B c +\frac {3}{8} x^{8} B a \,b^{2}+\frac {3}{10} x^{10} A a \,c^{2}+\frac {3}{10} x^{10} A \,b^{2} c +\frac {3}{5} x^{10} B a b c +\frac {1}{10} x^{10} B \,b^{3}+\frac {1}{4} x^{12} A b \,c^{2}+\frac {1}{4} x^{12} B a \,c^{2}+\frac {1}{4} x^{12} B \,b^{2} c +\frac {1}{14} x^{14} A \,c^{3}+\frac {3}{14} x^{14} B b \,c^{2}+\frac {1}{16} B \,c^{3} x^{16}\) \(194\)
parallelrisch \(\frac {1}{2} a^{3} A \,x^{2}+\frac {3}{4} x^{4} A \,a^{2} b +\frac {1}{4} x^{4} B \,a^{3}+\frac {1}{2} x^{6} A c \,a^{2}+\frac {1}{2} x^{6} A a \,b^{2}+\frac {1}{2} x^{6} B \,a^{2} b +\frac {3}{4} x^{8} A a b c +\frac {1}{8} x^{8} A \,b^{3}+\frac {3}{8} x^{8} a^{2} B c +\frac {3}{8} x^{8} B a \,b^{2}+\frac {3}{10} x^{10} A a \,c^{2}+\frac {3}{10} x^{10} A \,b^{2} c +\frac {3}{5} x^{10} B a b c +\frac {1}{10} x^{10} B \,b^{3}+\frac {1}{4} x^{12} A b \,c^{2}+\frac {1}{4} x^{12} B a \,c^{2}+\frac {1}{4} x^{12} B \,b^{2} c +\frac {1}{14} x^{14} A \,c^{3}+\frac {3}{14} x^{14} B b \,c^{2}+\frac {1}{16} B \,c^{3} x^{16}\) \(194\)
default \(\frac {B \,c^{3} x^{16}}{16}+\frac {\left (A \,c^{3}+3 B b \,c^{2}\right ) x^{14}}{14}+\frac {\left (3 A b \,c^{2}+B \left (a \,c^{2}+2 b^{2} c +c \left (2 a c +b^{2}\right )\right )\right ) x^{12}}{12}+\frac {\left (A \left (a \,c^{2}+2 b^{2} c +c \left (2 a c +b^{2}\right )\right )+B \left (4 a b c +b \left (2 a c +b^{2}\right )\right )\right ) x^{10}}{10}+\frac {\left (A \left (4 a b c +b \left (2 a c +b^{2}\right )\right )+B \left (a \left (2 a c +b^{2}\right )+2 b^{2} a +c \,a^{2}\right )\right ) x^{8}}{8}+\frac {\left (A \left (a \left (2 a c +b^{2}\right )+2 b^{2} a +c \,a^{2}\right )+3 B \,a^{2} b \right ) x^{6}}{6}+\frac {\left (3 A \,a^{2} b +B \,a^{3}\right ) x^{4}}{4}+\frac {a^{3} A \,x^{2}}{2}\) \(226\)

input
int(x*(B*x^2+A)*(c*x^4+b*x^2+a)^3,x,method=_RETURNVERBOSE)
 
output
1/2*a^3*A*x^2+(3/4*A*a^2*b+1/4*B*a^3)*x^4+(1/2*A*c*a^2+1/2*A*a*b^2+1/2*B*a 
^2*b)*x^6+(3/4*A*a*b*c+1/8*A*b^3+3/8*a^2*B*c+3/8*B*a*b^2)*x^8+(3/10*A*a*c^ 
2+3/10*A*b^2*c+3/5*B*a*b*c+1/10*B*b^3)*x^10+(1/4*A*b*c^2+1/4*B*a*c^2+1/4*B 
*b^2*c)*x^12+(1/14*A*c^3+3/14*B*b*c^2)*x^14+1/16*B*c^3*x^16
 
3.1.97.5 Fricas [A] (verification not implemented)

Time = 0.26 (sec) , antiderivative size = 166, normalized size of antiderivative = 1.00 \[ \int x \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{16} \, B c^{3} x^{16} + \frac {1}{14} \, {\left (3 \, B b c^{2} + A c^{3}\right )} x^{14} + \frac {1}{4} \, {\left (B b^{2} c + {\left (B a + A b\right )} c^{2}\right )} x^{12} + \frac {1}{10} \, {\left (B b^{3} + 3 \, A a c^{2} + 3 \, {\left (2 \, B a b + A b^{2}\right )} c\right )} x^{10} + \frac {1}{8} \, {\left (3 \, B a b^{2} + A b^{3} + 3 \, {\left (B a^{2} + 2 \, A a b\right )} c\right )} x^{8} + \frac {1}{2} \, {\left (B a^{2} b + A a b^{2} + A a^{2} c\right )} x^{6} + \frac {1}{2} \, A a^{3} x^{2} + \frac {1}{4} \, {\left (B a^{3} + 3 \, A a^{2} b\right )} x^{4} \]

input
integrate(x*(B*x^2+A)*(c*x^4+b*x^2+a)^3,x, algorithm="fricas")
 
output
1/16*B*c^3*x^16 + 1/14*(3*B*b*c^2 + A*c^3)*x^14 + 1/4*(B*b^2*c + (B*a + A* 
b)*c^2)*x^12 + 1/10*(B*b^3 + 3*A*a*c^2 + 3*(2*B*a*b + A*b^2)*c)*x^10 + 1/8 
*(3*B*a*b^2 + A*b^3 + 3*(B*a^2 + 2*A*a*b)*c)*x^8 + 1/2*(B*a^2*b + A*a*b^2 
+ A*a^2*c)*x^6 + 1/2*A*a^3*x^2 + 1/4*(B*a^3 + 3*A*a^2*b)*x^4
 
3.1.97.6 Sympy [A] (verification not implemented)

Time = 0.04 (sec) , antiderivative size = 199, normalized size of antiderivative = 1.20 \[ \int x \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {A a^{3} x^{2}}{2} + \frac {B c^{3} x^{16}}{16} + x^{14} \left (\frac {A c^{3}}{14} + \frac {3 B b c^{2}}{14}\right ) + x^{12} \left (\frac {A b c^{2}}{4} + \frac {B a c^{2}}{4} + \frac {B b^{2} c}{4}\right ) + x^{10} \cdot \left (\frac {3 A a c^{2}}{10} + \frac {3 A b^{2} c}{10} + \frac {3 B a b c}{5} + \frac {B b^{3}}{10}\right ) + x^{8} \cdot \left (\frac {3 A a b c}{4} + \frac {A b^{3}}{8} + \frac {3 B a^{2} c}{8} + \frac {3 B a b^{2}}{8}\right ) + x^{6} \left (\frac {A a^{2} c}{2} + \frac {A a b^{2}}{2} + \frac {B a^{2} b}{2}\right ) + x^{4} \cdot \left (\frac {3 A a^{2} b}{4} + \frac {B a^{3}}{4}\right ) \]

input
integrate(x*(B*x**2+A)*(c*x**4+b*x**2+a)**3,x)
 
output
A*a**3*x**2/2 + B*c**3*x**16/16 + x**14*(A*c**3/14 + 3*B*b*c**2/14) + x**1 
2*(A*b*c**2/4 + B*a*c**2/4 + B*b**2*c/4) + x**10*(3*A*a*c**2/10 + 3*A*b**2 
*c/10 + 3*B*a*b*c/5 + B*b**3/10) + x**8*(3*A*a*b*c/4 + A*b**3/8 + 3*B*a**2 
*c/8 + 3*B*a*b**2/8) + x**6*(A*a**2*c/2 + A*a*b**2/2 + B*a**2*b/2) + x**4* 
(3*A*a**2*b/4 + B*a**3/4)
 
3.1.97.7 Maxima [A] (verification not implemented)

Time = 0.19 (sec) , antiderivative size = 166, normalized size of antiderivative = 1.00 \[ \int x \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{16} \, B c^{3} x^{16} + \frac {1}{14} \, {\left (3 \, B b c^{2} + A c^{3}\right )} x^{14} + \frac {1}{4} \, {\left (B b^{2} c + {\left (B a + A b\right )} c^{2}\right )} x^{12} + \frac {1}{10} \, {\left (B b^{3} + 3 \, A a c^{2} + 3 \, {\left (2 \, B a b + A b^{2}\right )} c\right )} x^{10} + \frac {1}{8} \, {\left (3 \, B a b^{2} + A b^{3} + 3 \, {\left (B a^{2} + 2 \, A a b\right )} c\right )} x^{8} + \frac {1}{2} \, {\left (B a^{2} b + A a b^{2} + A a^{2} c\right )} x^{6} + \frac {1}{2} \, A a^{3} x^{2} + \frac {1}{4} \, {\left (B a^{3} + 3 \, A a^{2} b\right )} x^{4} \]

input
integrate(x*(B*x^2+A)*(c*x^4+b*x^2+a)^3,x, algorithm="maxima")
 
output
1/16*B*c^3*x^16 + 1/14*(3*B*b*c^2 + A*c^3)*x^14 + 1/4*(B*b^2*c + (B*a + A* 
b)*c^2)*x^12 + 1/10*(B*b^3 + 3*A*a*c^2 + 3*(2*B*a*b + A*b^2)*c)*x^10 + 1/8 
*(3*B*a*b^2 + A*b^3 + 3*(B*a^2 + 2*A*a*b)*c)*x^8 + 1/2*(B*a^2*b + A*a*b^2 
+ A*a^2*c)*x^6 + 1/2*A*a^3*x^2 + 1/4*(B*a^3 + 3*A*a^2*b)*x^4
 
3.1.97.8 Giac [A] (verification not implemented)

Time = 0.32 (sec) , antiderivative size = 193, normalized size of antiderivative = 1.16 \[ \int x \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=\frac {1}{16} \, B c^{3} x^{16} + \frac {3}{14} \, B b c^{2} x^{14} + \frac {1}{14} \, A c^{3} x^{14} + \frac {1}{4} \, B b^{2} c x^{12} + \frac {1}{4} \, B a c^{2} x^{12} + \frac {1}{4} \, A b c^{2} x^{12} + \frac {1}{10} \, B b^{3} x^{10} + \frac {3}{5} \, B a b c x^{10} + \frac {3}{10} \, A b^{2} c x^{10} + \frac {3}{10} \, A a c^{2} x^{10} + \frac {3}{8} \, B a b^{2} x^{8} + \frac {1}{8} \, A b^{3} x^{8} + \frac {3}{8} \, B a^{2} c x^{8} + \frac {3}{4} \, A a b c x^{8} + \frac {1}{2} \, B a^{2} b x^{6} + \frac {1}{2} \, A a b^{2} x^{6} + \frac {1}{2} \, A a^{2} c x^{6} + \frac {1}{4} \, B a^{3} x^{4} + \frac {3}{4} \, A a^{2} b x^{4} + \frac {1}{2} \, A a^{3} x^{2} \]

input
integrate(x*(B*x^2+A)*(c*x^4+b*x^2+a)^3,x, algorithm="giac")
 
output
1/16*B*c^3*x^16 + 3/14*B*b*c^2*x^14 + 1/14*A*c^3*x^14 + 1/4*B*b^2*c*x^12 + 
 1/4*B*a*c^2*x^12 + 1/4*A*b*c^2*x^12 + 1/10*B*b^3*x^10 + 3/5*B*a*b*c*x^10 
+ 3/10*A*b^2*c*x^10 + 3/10*A*a*c^2*x^10 + 3/8*B*a*b^2*x^8 + 1/8*A*b^3*x^8 
+ 3/8*B*a^2*c*x^8 + 3/4*A*a*b*c*x^8 + 1/2*B*a^2*b*x^6 + 1/2*A*a*b^2*x^6 + 
1/2*A*a^2*c*x^6 + 1/4*B*a^3*x^4 + 3/4*A*a^2*b*x^4 + 1/2*A*a^3*x^2
 
3.1.97.9 Mupad [B] (verification not implemented)

Time = 0.03 (sec) , antiderivative size = 169, normalized size of antiderivative = 1.02 \[ \int x \left (A+B x^2\right ) \left (a+b x^2+c x^4\right )^3 \, dx=x^8\,\left (\frac {3\,B\,c\,a^2}{8}+\frac {3\,B\,a\,b^2}{8}+\frac {3\,A\,c\,a\,b}{4}+\frac {A\,b^3}{8}\right )+x^{10}\,\left (\frac {B\,b^3}{10}+\frac {3\,A\,b^2\,c}{10}+\frac {3\,B\,a\,b\,c}{5}+\frac {3\,A\,a\,c^2}{10}\right )+x^4\,\left (\frac {B\,a^3}{4}+\frac {3\,A\,b\,a^2}{4}\right )+x^{14}\,\left (\frac {A\,c^3}{14}+\frac {3\,B\,b\,c^2}{14}\right )+x^6\,\left (\frac {B\,a^2\,b}{2}+\frac {A\,c\,a^2}{2}+\frac {A\,a\,b^2}{2}\right )+x^{12}\,\left (\frac {B\,b^2\,c}{4}+\frac {A\,b\,c^2}{4}+\frac {B\,a\,c^2}{4}\right )+\frac {A\,a^3\,x^2}{2}+\frac {B\,c^3\,x^{16}}{16} \]

input
int(x*(A + B*x^2)*(a + b*x^2 + c*x^4)^3,x)
 
output
x^8*((A*b^3)/8 + (3*B*a*b^2)/8 + (3*B*a^2*c)/8 + (3*A*a*b*c)/4) + x^10*((B 
*b^3)/10 + (3*A*a*c^2)/10 + (3*A*b^2*c)/10 + (3*B*a*b*c)/5) + x^4*((B*a^3) 
/4 + (3*A*a^2*b)/4) + x^14*((A*c^3)/14 + (3*B*b*c^2)/14) + x^6*((A*a*b^2)/ 
2 + (A*a^2*c)/2 + (B*a^2*b)/2) + x^12*((A*b*c^2)/4 + (B*a*c^2)/4 + (B*b^2* 
c)/4) + (A*a^3*x^2)/2 + (B*c^3*x^16)/16